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Abstract 

We study radion phenomenology in the context of flavor shining in warped extra dimension mod- 
els. In this unique setup, originally proposed by Rattazzi and Zaffaroni, solutions to the gauge hierar- 
chy problem and the new physics flavor problem are unified. A special role is played by the vacuum 
energy on the branes, that naturally allows for flavon stabilization and parametrically raises the ra- 
dion mass. We note that the radion mass squared is suppressed only by the log of the weak-Planck 
hierarchy, and it is in the favored range of the standard model Higgs. We emphasize that the radion 
to di-photon, to rr and to WW* can be promising discovery channels at the LHC, with a rate above 
that of the standard model Higgs. We find that the radion is unlikely to account for the excess in W 
plus dijet events as recently reported by the CDF collaboration. 

1 Introduction 

The Randall-Sundrum (RS) framework consists of a slice of anti de Sitter space in five dimensions 
(AdSs), where the warped geometry naturally generates the weak-Planck hierarchy Q and offers new 
approaches to flavor physics [|2]. It may also address the standard model (SM) flavor puzzle via the split- 
fermion mechanism [31, using flavor dependent wave function localization for the SM fermions dH [51. 
In addition, a protection from large flavor and CP violation is obtained via the so called RS-GIM mecha- 
nism [i4l[5l[6l[71[8l, even if the fundamental flavor parameters are anarchic. However, a residual little CP 
problem, in the form of sizable contributions to ex, e' /^k [|9l[lQl[IIl[l2l[T3l[l4]| and to the electric dipole 
moments [|71 [8l HSl [161, ruins such a nice feature. For further discussions on the RS flavor problems see 
e.g. [|l7l[Il[IllMinil2l|2llMl25llMl22ll2l. This RS Httle CP problem pushes the Kaluza Klein 
(KK) scale to be larger than C(10 TeV) |fT6l , which is beyond of the LHC reach and implies a severe little 
hierarchy problem. Gauging the SM approximate flavor symmetries (or part of them) in the bulk may 
solve this problem as was recently investigated by various authors [fT7ll25l l30l[3Tl [32ll33l[34l[35ll. 

In a large class of 5D flavor models the flavor symmetry is broken by unspecified dynamics on the 
UV-brane, and flavor violation is shined [|36l |29l [37l [381 to the IR-brane by scalar flavon fields, which 
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acquire vacuum expectation values (VEVs). It is interesting that a model of this type, where the full SM 
flavor hierarchies are generated by some mechanism on the UV-brane, can ameliorate the little hierarchy 
problem and improve the visibility of the model [fT6l [38ll . The above possibility is also motivated by 
the AdS/CFT correspondence [|39ll40ll . where the RS setup can be understood from a 4D point of view 
[|29ll4Tll42ll43ll44l| . identifying 4D global currents with bulk gauge symmetries. Indeed, in the context 
of electroweak precision tests, a significant improved Electroweak fit is obtained when the custodial 
symmetry of the SM is gauged in the bulk [|45ll . 

The gauge hierarchy problem can be addressed by the RS framework only if the size of the extra 
dimension is stabilized to (9(40)/A;, where k is the AdSs curvature. Therefore, the radion field, which 
corresponds to fluctuations in the distance between the branes, should develop a VEV (and a positive 
mass squared) that sets the size of the extra dimension. The stabihzation problem in the RS framework 
was studied by||4B|471|48l|43|53|5Il[51|53IM[55l|5^ when the favored solution is 

the Goldberger-Wise (GW) mechanism [|46ll47ll . It requires the presence of an additional bulk scalar field 
with a non-trivial background along the extra dimension, and a relatively small 5D mass. It is interesting 
to unify the fields which induce shining of flavor violation and stabilization, such as the GW scalar and 
the flavons. The idea of utilizing flavons as GW scalars was raised by Rattazzi and Zaffaroni (29] for 
IR localized SM model. Here we consider the case of ll37l . where the SM fermions and gauge fields 
propagate in the bulk, that allows one to solve the little hierarchy problem, and show that stabilization 
can be naturally induced by the flavons. 

Radius stabilization is induced by the presence of flavon VEVs and brane tension terms. Further- 
more, the radion mass squared is parametrically enhanced, as it is suppressed only by the log of the 
weak-Planck hierarchy, in agreement with the recent result of [|56l . The radion mass is found to be in 
the favored range of the SM Higgs mass. Hence, in this model the radion is likely to be the lightest new 
particle and may be the first signal of this framework. 

The remaining of the paper is organized as follows. In section|2]we study the GW mechanism in the 
presence of brane tension deviation and estimate the radion mass, while in section|3]we utilize the flavon 
fields as GW scalars. Section |4] contains a discussion on the collider phenomenology of the radion. We 
end with our conclusions in section |5]. The Appendices contain technical derivations. 

2 The Goldberger-Wise Mechanism and the Radion Mass 

The GW mechanism provides a classical solution to the stabilization problem of the RS framework. In 
order to ensure that the main appealing features of the RS framework are retained, we consider a small 
back-reaction of the GW scalar on the metric. The background solution of the gravity-scalar system is 
obtained using perturbation theory, following the method proposed in [fSTll . For simplicity, we use the 
4D effective potential method [|46l l47l . where the extra dimension is integrated out. In Appendix |C[ it 
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is shown that the result for the radion mass in this method is identical to the one that is obtained in the 
method of Il50l|58l. 



2.1 The Goldberger-Wise Setup 

We focus on the original RS setup and the GW mechanism. The extra dimension, y E {—yc,yc], is 
bounded by two 3 + 1 branes. One brane at y = and is called the UV-brane, and the other one at y = yc 
is called the IR-brane. Since we look for a background solution that takes into account the back-reaction 
of the GW scalars on the metric, and preserves 4D Poincare invariance, we choose the following ansatz 
for the metric 

ds^ = GMNdx^dx^ = Q-^^^y\^^dx^'dx'' - dy^ . (1) 

Greek indices run over (0, 1, 2, 3), capital Latin for (0, 1, 2, 3, 4), ?7^j, = diag(l, —1, —1, —1) and xa = y- 
We identify y with —y via orbifold, S1/Z2, assignment. The gravity part of the action forces an AdSs 
background metric, with curvature k. It includes the Ricci scalar TZ, a negative bulk cosmological con- 
stant Abuik = — 12M^/c^ and the brane tensions, which are parameterized by the dimensionless param- 
eters TiR/uv- The 5D Planck scale is M ~ 10^^ GeV and we assume that k < 2.6M (see discussion 
on section |4]). The GW mechanism contains bulk scalars, (pi, with bulk and brane potentials, which are 
^buik ^^'^ ^iR/uv(*^^) respectively, i will become flavor indices later on when we consider flavon 
fields as GW scalars. The considered action is 



^ ~ I ^^-^ I -^gravity + ^ ] -^GW ) ' 



where 



gravity = ^ {UM^^ - M'^TZ) - ^^M^kT^jy5{y) - ^M^kTi^,5{y - y,) , 



g is the determinant of the metric and S'ir/uv are the determinants of the induced metric on the branes. 
The GW scalars potentials are 



= \hK,ik (0. - ri/iR,M3/2) , (3) 
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where ej, feiR/uv.i, t" and ^'iR/uv,j are dimensionless parameters, (6iR/uv,i) parametrize the bulk 
(branes) masses of the GW scalars in units of the curvature and we consider as small parameters. For 
simplicity, we consider non-interacting quadratic potentials for the GW scalars. Higher terms of the GW 
scalars are negligible in the perturbative expansion, for further discussion see [|29ll . The boundary con- 
ditions of the GW scalars are determined by the brane potentials, and can be Neumann for 6iR/uv,j = 
or modified Dirichlet for &iR/uv,i ^ oo, where the boundaries' value of the GW scalars are set by r and 
^iR/uv,i- The result for finite values of 6iR/uv,i interpolates between the above extremes. Our parame- 
terization for branes' potentials is chosen due to the following reason. For large &iR/uv,i the GW scalars 
VEVs on the branes are 0i(O) ~ M^/^rz/uv,i and 0i(yc) ~ M^^'^rui^^^i, which mean that the GW scalars 
branes' potentials contribution to the brane tensions is small in this case. As explained in section 2.2 , r 
serves as the expansion parameter for the background solution. Therefore, for later use we parametrize 
the brane tension as follows 

TiR = t/°) + r^^TiR , Tuv = T^'i + r^ST^y . (4) 

2.2 The Background Solution 

The background solution preserves 4D-Poincare invariance, i.e. we consider solutions where the GW 
scalars and the metric backgrounds are only y-dependent. The metric is symmetric under the orbifold 
symmetry and we assume that for the GW scalars as well. The GW scalars equation of motion and the 
5D Einstein equations, which are derived from Eqs. ([T])-(|3]) are 



Ky - Vc) , (5b) 



+ 



M^kTi^ + Ui^,k (0, - rz/iR,M3/2)^ 



^'=^' + ^A'i>?-^m , (5c) 



where repeated indices are summed over and ' are derivatives with respect to y. From Eq. ( pq ) one can 
see that non zero VEVs of the GW scalars will perturb the metric away from the pure AdSs form. Later 
we will see that the background solution of the GW scalars is proportional to r (and higher powers of it), 
therefore it naturally serves as an expansion parameter and we can ensure that the back-reaction of the 
metric from the GW scalars is small. 

The zeroth order equations can be obtained by setting r = and 4>i{y) = 0. After we integrate over 
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Eqs. ([5b]) and ([5cJ) (where Eq. is plugged into Eq. ([5b])) we get 



A{y) = k\y\, Ti°^ = 12, t/^^ = -12 . (6) 

This is the pure AdSs solution, and the values of Tj^^/uv those of the unstabilized RS model. There- 
fore, 5TiR/uv the brane tensions deviations. From now on, we set v4(0) = 0, which is a gauge fixing. 
In Eq. ([6]) we chose the positive solution of the square root, A = \k\ (the negative solution. A' = —\k\, 
corresponds to interchanging the IR and UV branes). 



To first order in r, Eq. ( [5a| ) in the bulk is 

(j)'! - Ak(f)', - eik^cj), = , (7) 

with the boundary conditions 

(0) = 6uv,^fc (<^i(0) - rz/uviM=^/2) , (8a) 
2cj)[{y,) = -hi^^.k {Uy,) - riyjn,,M^/^) . (8b) 

The solution of Eqs. is 

Uy) =rM^I^ (^Xi,e('-^)'=l^l + X2,e(2+^)'l^l) , (9) 
where Xj i and Xj 2 are given by 

Xii= ^ — , — , r , , (10a) 

e2fc,.V4+^7 (2^4+^ + 4 + 6ir) + 2v/4+^ - 61R - 4 

Xi^2 =i^uv,i - Xi^i ■ (10b) 

and consider the the limit of 6uv,j — ^ 00. As long as 6uv,j ^ taking finite 6uv,i does not change our 
results. As we claimed before, 0i oc r therefore A{y) and Tuv/ir receive no contributions at first order 
in r. 

To order the GW scalars' equations are left unchanged, but the one involving A{y) is 

^iy) = k+^^my)-ek^m) , (11) 

and after integration we get 

A{y) = k\y\+r^kG{\y\) , (12) 
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where 



- e^'^^v^Xi.iXa.e, + Xj, (2 - ^4+^) 



(13) 



G(?/) is the leading order correction to the AdSs metric from the GW scalars back-reaction. 

The expansion in power series of r naively requires r ^ 1, but in the present case one can relax this 
condition. In fact, for a sizable range of parameter space, discussed below, r can be of order unity. From 



Eq. (5c ) we find that 



A' = k\ 1 + 



0f - 



24M3A;2 



(14) 



In order to have an AdSs-like solution we expand the square root in the above equation to lowest non- 
trivial order. The condition for the validity of this approximation is 



L'2 



2 S/fS 



(15) 



which, as anticipated, allows for r ~ 0{1) for a wide range of the parameters. 



2.3 The 4D Effective Action and the Radion Mass 

The 4D effective action is derived by integrating out the fifth dimension, similarly to [|46l , for detailed 
derivation see Appendix [A|. Let us denote the radion field as lpr{x^) = fa{x^), where 



a[x' 



f 



12- 



M3 



(16) 



For successful stabilization the radion VEV should produce the weak-Planck hierarchy, therefore (a) = 
~ 10~^^. The effective action to order and leading order in e and a is 



-Ceff - Kfr(a) + const 



(17) 



where 



4b 



+ 2a'^'^uly^i (foi^R , - 32) - 4a^>/Vuv,i'^iR,i&iR,i(4 + 6ir,»; 



(^IR.i'^IR.i 



(18) 
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The constant in Eq. ( [17] ), sets the value of the potential to zero at its minimum, by tuning the value 
of (5Tuv- The corresponding fine-tuning in this procedure is nothing but the celebrated cosmological 



constant problem, whose solution is beyond the scope of this work. Eq. ( [18] ) gives the leading order form 
of the effective potential, we note that G{y) contributes only to order of the effective potential. 

The radion mass is derived by using the effective potential. For simplicity we consider the case of 
one GW scalar, and denote Cj = e, 6iR,i = 6ir, z/uv,j = 1 and z/jR^j = u. To leading order in e and 
1/ log(a^), the extremum condition, dV/da\ = 0, is 



a4 =u± 



-5Tir6ir(8 + 6ir) 



2b 



+ 0(e) 



IR 



which are the two extrema to of the effective potential. The radion mass is 



m 



1 d^V{a) 



1 



,^, 2 2 2 V-5Tir6ir(8 + 6ir) 



8 + 6 



-5Tir6ir(8 + 6ir) 



IR 



2b 



IR 



(19) 



(20) 



where the ± refer to the different solutions of -t respectively. From Eq. ( [20[ ) we can see that one 
extremum corresponds to a minimum of the potential and the other to a maximum. For e > and 
(5TiR < the minimum is at a^,^ for 6ir > and at for 6ir < —8. For (5Tir > and —8 < 6ir < 
the minimum is located at a*+. For e < the locations of the maximum and the minimum are reversed. 
The desired hierarchy can be achieved without fine-tuning of u, 6ir and 5Tir. 

Unlike previous works but in agreement with [[56l . we find that the parametric dependence of the 



radion mass is m^^^ ~ fc^a^e, and not m^^^ ~ fc^a^e^/^ as in [|47l [^or m'^^ ~ k'^ale^ as in fSO^I, where 
the IR-brane tension is tuned. This parametric enhancement is due to the presence of IR-brane tension 
terms, and this is a general consequence of the GW mechanism. The original result of [47], where 
'^rad ~ fc^a^e'^/^, can be easily reproduced in the limit that the IR-brane tension is fine tuned (|5TiR|^e). 
From Eq. ( |T9| ) one can see that e ~ l/log(a^); thus, the radion mass square is inversely suppressed by 
log of the weak- Planck hierarchy (log(a^)). For typical values of parameters we get that 



m^ad = C(lOOGeV) 



(21) 



Furthermore, our result is not affected by the Casimir energy which is induced by bulk fields, see Ap- 
pendix [Bj 

We can get some intuition for the radion mass estimation by invoking the AdS/CFT correspondence. 

'in BtII the radion mass squared that is given in the main text is m^^^ ^ k'^ale^. However as the authors mentioned in 
footnote 2, a careful calculation yields m^^^ ^ k'^a^e^^^. 
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The conformal field theory (CFT) is defined below some UV cutoff, ~ k, and is broken spontaneously 
at some low energy scale, /a^ ~ fca* [|29l HTll . The radion corresponds to a 4D dilaton, the Goldstone 
boson of the broken conformal symmetry, and each 5D bulk field has a dual operator in the 4D CFT The 
4D dual of the GW scalar is a scalar operator with scale dimension of 4 + e/4 for |e| <^ 1. This operator 
explicitly breaks the conformal symmetry. Therefore conformal invariance is only an approximation 
and the dilaton becomes a pseudo-Goldstone boson, with a small mass. The same happens in the 5D 
theory, where the radion acquires a mass due to the GW mechanism. The mass square of the dilaton 
is proportional to the IR breaking scale times the parameter that explicitly breaks the symmetry. This 
parameter is proportional to the GW scalar mass and is given by e/4. Therefore, naively the dilaton mass 
is m^ji ~ k'^ale, which agrees with the result in Eq. ( |20l ) up to order one coefficients. The same scaling 
of the dilaton mass can be also obtained by trace anomaly matching [61 J . 



2.4 Stabilization and Fine-Tunings 

As discussed in the introduction, the gauge hierarchy problem can be addressed by RS only in the 
presence of stabilization mechanism, which reduces the number of required fine-tunings from two in 
the unstabilized RS to one. The remaining one corresponds to the 4D cosmological constant prob- 
lem [|46ll57ll59ll . which is left unsolved. In the absence of stabilization, the radion 4D-effective potential 
is 

r3rr I ^bulk , A ( ii^j^rp ^bulk 



Kfr(a) = fcM^Tuv + ^ + [kM'Ti^ - -f^ j . (22) 

The two fine-tunings correspond to (i) having a flat potential for the radion, otherwise the radion's VEV 
wiU be driven to zero or to infinity depending on the sign of the term; and to (ii) the vanishing of the 
4D effective cosmological constant, Kff (a*) = 0. The resulting conditions are 

rri rriiO) Afaulk rri rTi{0) Afoulk „ 

^IR - ^IR - - -1^ . ^UV - JuV - - 1^ • (^^) 

However, a stabilization mechanism results in generating new terms in the effective potential. These lead 
to stable configuration with non trivial VEV. Consequently, the need for the first fine-tuning is eliminated. 
The second fine-tuning, that corresponds to the vanishing of the 4D-effective cosmological constant, is 
still unavoidable. 

The same is also manifested in the 5D picture, where the warp factor has to satisfy the jump condi- 
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tions. They are derived by integrating of Eq. ( [Sb] ) over a small interval around the branes, 



G'(0) = ^ . (24b) 



It is straightforward to show that the value of yc, which is obtained by the IR-jump condition in Eq. (24a), 



is equivalent to the value of yc that minimizes the radion 4D-effective potential and that the UV-jump 



condition from Eq. ( |24b[ ) is equivalent to the demand of the potential vanishing at its minimum. In the 
absence of the GW mechanism = G{y) = 0), the two jump conditions do not depend on y^ there- 
fore 5TiR/uv = is required. This corresponds to both of the above fine-tunings in the RS-framework. On 
the other hand, when adding the GW mechanism, the jump conditions depend on yc. Thus, manipulation 
of the value of yc to fulfill the IR-jump condition is possible, which allows also for brane tension devia- 
tion. Once yc is set, there is no such freedom in the UV-jump condition which is satisfied by fine-tuning 
of the UV-brane tension, again resulting with a single fine-tuning at the end Il57l . 



3 Flavon Stabilization, Flavons as Goldberger-Wise Scalars 

In this section we focus on the model of [|371 and concentrate on the quark sector only. In that model, 
the SM fermions and gauge fields propagate in the 5D bulk, while the SM Higgs is IR localized. An 
SU{3)q X SU{3)d subgroup of the SM flavor symmetry is gauged in the bulk and it is broken by bi- 
fundamental flavon field, y^. Although we focus on a specific model, our result also applies to other RS 
shining models. For example, in the model of [|38l , one can consider, for instance, the case of a small 
bottom Yukawa, where the main contribution to stabilization comes from the top flavon sector (then the 
bulk field is ?/„ corresponding to up-type Yukawa and not yd). In addition, our analysis can be extended 
to other models where both up and down Yukawa propagate in the bulk, such as ll38l with sizable bottom 
Yukawa, where the main contributions would be from the third generation sector. In these cases one 
expects an even richer potential and stabilization options to arise. 

The quantum numbers of the 5D counterparts of the SM doublets and down-type singlet quarks under 
the flavor group are \&q ~ (3, 1) and ~ (1, 3), while the flavon quantum numbers are yd ~ (3, 3). The 
5D fermions bulk masses, cx (in units of curvature), and the dimensionless Yukawa, Yd, axe to lowest 
order in the flavons VEV, (yd) 

CQ = aQl + PQ{yd){yl) = «q1 + PQYdYj , (25a) 
Cd = adl + My\) iVd) = «dl + /3dYjYd , (25b) 
kXd{yd) = Yd. (25c) 
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We assume that the flavor symmetry is not broken explicitly in the bulk. Therefore, to lowest order, the 
bulk action for the flavon is 



(26) 



where e is dimensionless and parametrizes the bulk mass of the flavon fields and we assume e ^ 1 as in 
section [2j It also ensures that the RS flavor solution is not ruined, as for small e the flavon profile along 
the extra dimension will be mostly flat, see Eq. (|9]). It is straightforward to show that for small e the flavor 
structure of bulk RS models retains its qualities even with non-constant bulk fermion masses resulting 
from the flavon fields [|37l . 

As discussed, the flavor symmetry is broken only on the UV-brane, such that the UV boundary 
condition for the equation of motion is 



^(diag) 



yd 



( 



\ 








\ 

^'UV,2 

i^UV,3 ) 







(27) 



We have parametrized the fields' values through the dimensionless parameters z^uv,*. which one can 
normalize by defining z^uv,3 = 1: and define i>uv = ('^uv,!)^ + ('^uv,2)^ + {y\^v£f ■ We assume that 
all of the entries are real. For the anarchic case, all the entries are of the same order of magnitude 
and therefore z>uv — V^, while for hierarchical cases the main contribution will come from the third 
generation entry, and thus z/uv ~ 1- Given the naive dimensional analysis estimates for the values of the 
flavons from [|37ll . and taking into account the factor two tuning of the {yd) found there, we find 



(28) 



where yU < M is the scale at which the theory becomes strongly coupled. Therefore, r < 0.4 and the 



perturbative expansion, developed in section 2.2, is expected to be valid. The UV boundary conditions 
can be formally written using the 6uv — ^ cxo limit of the UV brane potential 



\yd 



_ £)(diag) I 



(29) 



We assume that on the IR-brane the flavor symmetry is not broken explicitly. Therefore, the lowest order 
term which one could write is a brane mass term 



ydy\ 



(30) 



This is the only relevant term one could write for the brane potential, invariant under the flavor symmetry. 
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Once the potentials are specified, the background solution of the gravity- scalars system can be cal- 
culated perturbatively in r as in section 2.2 . The complex fields y'J are separated to real and imaginary 



components y^^ = {y^m + Wid) ■ the notations of section 2.2. the flavons' IR-brane potential corre- 
sponds to z/jR j = for all fields, while in the UV-brane potential z/uv,j 7^ for the real diagonal entries at 
the limit 6uv,j cxd (z^uv.i = for all the other fields). The resulting background solution vanishes for 
fields with zero UV boundary condition, which are the off-diagonal fields and imaginary diagonal fields, 
Vmiy) = for i 7^ j and yf^iv) = for V (i, j). For the real diagonal fields, the solution is 

VUV) = M'/\u^y, (^Xie^^-v^)'^!^! + X^e^^+^^^l^l) , (31) 

^ ^ e^^^^v^ (4 + 2v/4T^ + 6ir) X =1-X 

' e^f^y^^' (4 + 2v/4T^ + 6ir) - 4 + 2^/IT~e - 6ir ' ' ' ' 

which is the result given by Eq. (|9]) in the limit of i^ir j — > . The radion effective potential is calculated 
by integrating out the extra dimension. Since the three fields have the same profile but it is rescaled 
by different values on the UV-brane, and because of the fact that the effective potential is the result 



of integrating only terms quadratic in the fields, the effective potential is simply Eq. ( [T8| ), by setting 
J^iR 1 = 0. The result is 



2„4 



VeM =M'kr'a 



+ 0(e,a^), (32) 



where corrections from back-reaction of order have been neglected. 

In contrast to the general case of GW mechanism, this potential has only one extremum. It is 
important to note that fiavon stabilization may occur only with an IR-brane tension deviation term (notice 
that these results can change when sub-leading irrelevant operators on the IR brane are considered). The 
potential extremum, which is at a = a*, and the radion mass to leading order in e and are given by 



el2 '^^IR (8 + ^Ir) , r^i \ 2 1 r 2 2 2 c-TH , /o^ 2\ /oon 

aj =-—^ 77 + m^^^ = --k a^er 5Ti^ + 0[e ) , (33) 

z^uv 46iR 6 



where we assume 6ir ^ e. Because al^"^ ~ for the desired hierarchy and small bulk mass, it may 



be inferred from Eq. ( [33] ) that as in the original GW mechanism, the desired hierarchy can be reached 
without fine-tuning of 6ir and (5Tir. The radion mass scales as in the single GW scalar case considered 
above, with typical values of order (9(100 GeV). 

Note that for 6ir — 0, which is equivalent to imposing Neumann boundary conditions on the IR- 
brane, the mass vanishes regardless of the sign of Careful re-derivation of the radion mass, to second 



^ This result is quite general, and in fact one can show, using the effective potential method, that for GW scalars with 
Dirichlet and Neumann boundary conditions on the UV and IR branes respectively, there cannot be stabilization with the 
desired hierarchy for all values of e (with no assumption on its size) even with the inclusion of the IR-brane tension deviation. 
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order in e, shows that in order to have a positive mass squared, one needs 6ir > e/2 or 6ir < — 8 + e/2 for 
e > and —8 — |e|/2 < 6ir < — |e|/2 for e < 0. We can claim, therefore, that a necessary condition for 
stabilization in the order we work in r is that there is a non- vanishing IR-brane potential for the flavons. 

There are further constraints on the parameter space from flavor physics. The profile of the flavon 
fields depends on the values of 6ir and e. For small positive e the fields' value decreases slightly in the 
bulk when moving away from the UV-brane, while it increases for negative e. Its value on the IR-brane, 
which determines the 4D fermion masses, is determined largely by the value of For 6ir ±oo, 
the fields' value on the IR brane vanishes, which is problematic for models in which the Higgs is an IR 
brane locahzed field. Conversely, for 6ir ~ —8 the fields' value on the IR-brane diverges. This signals 
a breakdown of the perturbative expansion, and we expect that higher order terms in the brane potentials 
will be important for this parameter choice. 

One should note that gauge symmetries in the bulk and their breaking on the boundaries have a 
dual description using the AdS/CFT correspondence. Gauge symmetries in the bulk of unsliced AdS 
correspond to global symmetries of the dual CFT Adding the UV and IR branes means that the bulk 
gauge fields can be decomposed to zero and higher KK modes. The higher KK modes are interpreted as 
composite states of the strongly coupled CFT. While the existence of massless modes of the gauge fields 
means that the symmetry is also weakly gauged in the CFT. If the symmetry is broken in the gravity 
side by boundary conditions the zero mode will pick up a mass. If the mass is of the order of the CFT's 
UV cutoff, the zero mode decouples and the symmetry is in fact a global symmetry on the 4D side; on 
the other hand, if the mass is of the order of the IR cutoff, the interpretation is that the CFT breaking 
also dynamically breaks the symmetry. Therefore, imposing the flavor symmetry and breaking it on the 
UV-brane, as we do, implies that the flavor symmetry of the SM is imposed as a global symmetry of the 
CFT. 

The flavon fields are then interpreted as dual to scalar operators on the conformal side which have 
non trivial quantum numbers under the flavor symmetry. If the flavons have small bulk masses, adding 
these operators explicitly breaks the conformal symmetry only by marginally relevant or irrelevant oper- 
ators. The pre-factors of the two exponents in the solution of the scalar equation of motion are related 
to the source of the CFT operator and its VEV respectively 04 11 . If an IR-brane potential for the flavon 
field is not added, one finds that the solution is nearly constant in the bulk and that the coefficient of 
the exponent proportional to the VEV is negligible, implying that inner CFT dynamics do not break the 
flavor symmetry. On the other hand, turning on only the brane mass terms, as we did, also turns on the 
coefficient of the exponent proportional to the VEV of the scalar, signaling that now the CFT breaks the 
flavor symmetry dynamically. However, because in this case the flavon VEVs do not twist in flavor space 
while propagating in the bulk [|62l . one interprets that the CFT dynamical breaking is aligned with the 
explicit breaking of the flavor symmetry, hence no additional flavor violation is induced. 
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4 Collider Phenomenology 



We focus on the case of light radion, which in the CFT description can be viewed as a pseudo Goldstone 
boson of the broken conformal symmetry. In this case, as discussed above, the typical value of the radion 
mass is (9(100 GeV), which interestingly coincides with the preferred range of the SM Higgs mass [|63l . 
We compare the radion signal to that of the SM Higgs in three interesting channels, gg ^ r / hgu 11, 
gg —7- r/hsM t^t^ and gg — )■ r/hgu WW*. We emphasize that in the first channel the radion 
signal can be naturally enhanced by a factor of (9(10). In the two other channels the radion signal can 
be larger by a factor of (9(5) compared to the SM Higgs case. In addition, we show that the radion is 
unlikely to account for the excess in the W plus dijet events as recently reported by the CDF collaboration 
[|64l[65||. It is important to note that in an equivalent study of the DO collaboration there is no excess [|66l . 

The radion collider phenomenology has been discussed in details by [|471 [50l EH |671 |68l |69l USi IZll 
1721 1731 174l 1751 1761 1771 ITU |79l- Here we focus some of the important features relevant to our analysis, 
especially in view of the LHC improved sensitivity to the above class of observables. The radion cou- 
plings to the matter fields can be obtained from general principles, similar to the known bonafide dilaton. 
The couplings are proportional to the mass of the fields (or more precisely to the effective 4D trace of the 
energy momentum tensor when integrated over the extra dimension [75]), 



Arad is the most important parameter for the radion phenomenology, and can be interpreted as a symmetry 
breaking scale in the 4D CFT language [1751 . The value of Arad is roughly around (9(TeV), but it 
can be viewed as a free parameter of the theory as long as it satisfies the perturbative bound [|80l [8TI . 
The perturbative limit implies that fc/Mpi < 3, where Mpi is the 4D Planck scale. From the relation, 
M|i ^ 2M^/k (the notation of ^ is used), we obtain that k/M < 2.6. Therefore, for ka^ ~ (9(1 TeV), 
there is no problem to consider A^ad ~ (9(1 TeV). 

In the (9(100 GeV) mass range, just as in the SM Higgs case, the di-photon discovery channel is 
one of the most important channels for the radion search at the LHC. In fact, it is found that the radion 
discovery potential in the (7(7 — r — )■ 77 process can be enhanced compared to gg — hsu -> 11 
for the reasonable theory parameter space [75]. The reason for the enhancement is as follows. The 
dominant contribution to the SM Higgs production via the gluon fusion involves triangle loop diagram 
with heavy fermions, and therefore, is proportional to the QCD beta function coefficient of the top quark 
(bt) [l82l . However, the dominant contribution to the gluon fusion production for the radion comes from 



(34) 



where i stands for a generic matter field, and Arad is defined as 




(35) 



13 



the trace anomaly, i.e., it is (in the CFT language) proportional to the beta function coefficient of massive 
composite fields (bcFx) [VS]. Because bcFT ^ bt, it follows that the radion coupling to a pair of gluons 
can be larger than that of the SM-Higgs. On the other hand, the branching ratio of the radion to di-photon 
is similar to that of the SM Higgs in most of the parameter space [|75l . 

In order to compare the di-photon channel of the radion and the SM Higgs at the colliders, it is 
useful to define an approximate formula for the ratio of the discovery significance of the radion in the 
^f^f — )• r — )• 77 channel and that of the SM-Higgs with the same mass, similar to the definition given 
in Bia, as 

P77 - ^(^) r(r -> gg) B{r 77) 



S{r) ^ r(r -> gg) B{r 77) 
S{hsM) r(/isM -> gg) B{hsM 77) ' 



A similar ratio can be defined for gg ^ r ^ t^t^ and gg ^ r WW* channels. Note that since both 
the SM Higgs and the radion are of narrow widths in the range of our interest, we do not include invariant 
mass resolution effects. By a careful examination of we can see that the ratio of the widths weakly 
depends on the radion and the SM Higgs masses and it has a simple Arad dependence 

r(r gg) r ( ^ B{r ^77) 



where v ~ 246 GeV is the SM Higgs VEV and /i(mrad) is monotonically changed by C(10%) in the 
interval of (9(100 GeV) which covers the mass range we are interested in. /2(mrad) is a non-trivial 
function of mrad and its value varies within 0{\). In Fig.[T](LHS) we show that the signal of the radion 
can be easily enhanced by (9(10) compared to the SM Higgs with the same mass. We assume that 
there is no Higgs-curvature coupling, and thus Higgs-radion mixing is absent. We also do not include 
brane-localized kinetic terms. In Fig.[T](RHS) we show the differential cross-section, da/dM^^, for the 
LHC (at center of mass energy of 7TeV) as a function of the di-photon invariant mass for a radion vs. 
SM Higgs of a mass of mrad,/isM — 120 GeV and Arad = ITeV . The cross section is simulated by 
MadGraph/MadEvent 5.1 [[831 [H and we use CTEQ6M PDF set. 

Another two interesting channels for the LHC are gg r ^ t^t^ and gg r ^ WW*, where 
the discovery significance is defined similar to the one of 77 in Eq. (|36]) and are labeled Rl^^'^ and R^^^* 
respectively. The general formula for the coupling of the radion to the SM-fermions with the bulk mass 
parameters cl and cr, defined in [|75l . is 



1 m/ (1 + 2cl) + al'^^+^'^jl - 2cr) + 2al+^'^^{cR - cl) - 2af^-^ 

2A,ad (1 - a^-^) (1 - ai+^^«) " ^''^ 

InFig.|2|(LHS) we show that the signal of the radion can be easily enhanced by (9(5) compared to that of 
the SM-Higgs for a given mass. As we can see the signal significance is very sensitive to the localization 
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Figure 1: Left: contours of ratio of discovery significance for a radion vs. SM Higgs of the same mass, 
FCg , in the plane of mrad-Arad- Right: da/dM^^ for a radion vs. SM-Higgs of the same mass for 
m^ad = 120 GeV and A^ad = 1 TeV. 

of both LH and RH r's in the bulk, i.e. cl and cr. In Fig.|2](RHS) we show that the signal of the radion 
to WW* can be bigger by a factor of 0{5) compared to the SM Higgs signal. 




-1.0 -0.5 0.0 0.5 1.0 130 135 140 145 150 155 
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Figure 2: Left: contours of ratio of discovery significance for a radion vs. SM Higgs of the same mass, 
R's^'^ in the cl-cr plane for rrij-ad = 120 GeV and Arad = 1 TeV. Right: contours of ratio of discovery 
significance for a radion vs. SM Higgs of the same mass, R^^\ in the plane of m^ad-Arad 



In a recent CDF study there is 4. la excess in the Wjj event sample in 7.3 fb^^ of data for dijet 
invariant masses between 120 — 160 GeV, which corresponds to a sizable cross section of 4pb [|64ll65l . 
However, no excess was found in a similar study done by the DO collaboration [|66ll . Although the radion 
mass can be naturally in the range to explain this anomaly, we find that it is unlikely to account for it. 
Below, we analyze the case where the SM fields propagate in the bulk. It is worth emphasizing that this 
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analysis holds also for the original RS model (where the SM is IR-brane localized). This is due to the 
fact that the correction to the ly+iy -radion coupling is about only 20% compared to the bulk SM case 
[|75l and therefore will not change our conclusion. 

The dominant channel for Wjj production that involves the radion is Higgs-like radion-strahlung. 
Other channels are proportional to light quark masses and are negligible. Since the phenomenology of the 
radion and the SM Higgs are similar, it is easy to estimate the ratio between the rate for pp — )■ Wr — )■ Wjj 
(dominanted by jj = gg) and pp — )■ Whsu — ^ Wjj (dominanted by jj = bb), which it turns out to be 



where ^ 1.4 comes from the sub-leading corrections to the iy+iy~-radion coupling, with v/Aj-ad ~ 
1/4 and BR{r — )■ jj)/ BR(hsM — ^ jj) ~ 1- Explicit MadGraph/MadEvent simulation leads to 
a{pp ^ Wr Wjj) ~ 1.1 fb, while a{pp WhsM)BR{h^M bb) ~ 12fb [85], which is con- 
sistent with our estimation. Therefore, we conclude that the Wjj anomaly can not be explained by the 
RS radion (contrary to the statement of [|86l ). 

5 Conclusions 

Flavon stabilization of warped extra dimension form a unique mechanism, originally proposed by Rat- 
tazzi and Zaffaroni, where solutions to the gauge hierarchy problem and the new physics flavor problem 
are unified. As we have shown in this setup, the stabilization is induced by a combination of flavon vac- 
uum expectation values and the presence of IR-brane tension term. Parametric enhancement of the radion 
mass is induced by the brane vacuum energy. It scales as ~ k^ale, in agreement with [|56l . where 
e ~ l/log(av,) ~ 1/40, corresponds to the weak-Planck hierarchy and ka^, ~ C(TeV). Therefore, 
the natural range for the radion mass is C(100 GeV) which is in the favored range of the standard model 
(SM) Higgs mass. Note that the radion mass is large enough to avoid radion-mediated flavor changing 
neutral currents [87], but still small enough to induce a sizable Sommerfeld enhancement that contributes 
to the dark matter (DM) annihilation, which may be probed via indirect astrophysical signals for TeV 
scale dark matter in Randall-Sundrum models ll88l |5 

In this mass range the di-photon channel, gg ^ r / hsu 11, forms an important discovery channel 
for both the radion and the SM Higgs, which tend to have similar collider phenomenology. We have 
emphasized that in this channel the radion signal can be naturally enhanced by C(10) compared with 
that of the Higgs. We also pointed out that in two other interesting channels, gg r/hsu ~^ r^r^ and 

^In RS framework, the first DM model was based on a grand unified theory (GUT) model 1891 190| . where stability of 
the DM is a spin-off of suppressing proton decay, and later extended to incorporate custodial Z ^ bb into the GUT frame 
work [88 1 . It is clear that radion mediated Sommerfeld enhancement might also be relevant for other RS-type scenarios with 
a DM candidate ll9Tll92ll93l . with a exception for 1941 . where DM candidate is KK-odd radion. 




(39) 
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gg — 7- r/hsM WW*, the radion signal is similarly enhanced by a factor of 0{5). We find that the 
radion is unlikely to account for anomalies in the W plus dijet differential distribution. 
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A Derivation of the Radion Effective action 

Here we derive the radion effective action following the method of ^46[ l47ll . The radion is identified 
with a simple field (x^ dependent) parameterization of the fifth dimension radius. The fifth dimension 
coordinate is transformed into polar coordinate, y = T{x^)0 = TO and E {—n,n], where the line 
element is given by 

ds^ = Q-^^^'^\'\\^,dx^'dx^ - T^dO^ . (40) 



and A{T \ 9\)h given by Eq. ([T2j). The 4D-effective action is derived by integrating out the fifth dimension 



^Sgw = f '^A ^g-ity m + C'S^ [0, {TO)] \ . (41) 



where 



^gravity = M^Tc"^^ ^12^ - 7^ - ^ [5(^)Tuv + 5{e - 7r)riR]| , 

7^Te-^^ = 6e-2A (^2y^/2 _ ^ ^/^ ^^^^^2 ^ ^^-iA ^^j.j^2 _ ^^a" _ 4^' {^^{6) - 6{e - vr))] . 

To order and leading order in a^, and e (note that we do not expand terms such as q-^'^'^'^ because 
k7re{T) ~ (9(1)), by using Eqs. ^ and ([16]) the 4D effective action is given by 

^Sgw =lf (1 + («)) id.af - Kfr(a) + const . (42) 
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where Ves{a) is given by Eq. ( [18] ) and 

Since we look for a solution around the VEV of a, denote (a) = ^ 1 which is independent, 
h{a) can not change the value of a^. Therefore, we plug a^, into /i(a)|^and omit all the interaction terms 
which are proportional to 9^a. Schematically, the 4D action is of the form 

^ts,GW = If fl + ^^Ha.)) {d,af - r'V{a) , (43) 



2'' V 12 



where r'^V{a) = Ves{a). Moving to a canonical basis where the appropriately normalized radion field is 
denoted by ipn, £efr,GW is given by 



2 \/y/l + p(a. 

= ^(5^Rf-r^^(y)+(^(r^)- (44) 

As one can see the contribution of h{a^ vanishes at order r^. The canonical form of the radion field is 



^K(a:'^) = /a=^^e-'=-^(^^). (45) 

B Goldberger-Wise and Casimir Energy Stabilizations 

A different approach from GW-mechasim to the radion stabilization is to implement the Casimir energy 
which is induced by bulk fields [|52l [531 |55l. Although the Casimir energy looks as a natural candidate 
for stabilization because one does not need to introduce a special field for stabilization (as in GW mech- 
anism), it is unnatural since it involves a small parameter, (9(10~^) or smaller, in order to get the right 
hierarchy Il53ll55l . This small parameter induces a fine-tuning, and therefore this mechanism can not be 
viewed as a complete solution to the hierarchy problem. However, in models that contains GW mecha- 
nism and bulk fields, one can ask if the Casimir energy ruins or reinforces the stabilization caused by the 
GW mechanism. 

In order to address this question we consider the GW- stabilization and the dominant contributions to 
the Casimir energjj^from bulk gauge fields. The Casimir energy due the bulk gauge fields were derived 

^For small e, the numerical value of hi^a^}) is of order unity. 

^BuUc scalars with mass = —Ak^ and bulk fermions with bulk mass, m/ = fc/2, induce similar contribution to the 
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in [|52l . In the case of one GW scalar and bulk gauge fields the radion 4D-effective potential is: 



VM=r'kM'aU6T,^ + 



46 



IR 



7^/4 _ 



4(8 + bm) 



8 + 

+ 2a^/2 - 32) - 4a^/V6iR(4 + bjn) 



1,2 2 



+ 



16n^ log(a) 



(46) 



where 



/3(Pir) 



dtf 



t^PmIi{t)+tIo{t) 



Ki{t) and Ii{t) are Bessel functions, Ng = QpN, gp is the number of the physical polarizations (gp = 3 for 
5D bulk gauge fields) and pi^ is the coefficient of the IR-brane kinetic term of the gauge fields. /3(pir) is 
an C(l) parameter and its sign depends on the value of piR (for example /3(0) = 1.005, /3(5) = —1.519 
and /3(oo) = —2.330). The extremum and the radion induced mass from the potential in Eq. ( |46l ), 
to leading order in e and 1/ log(a^), are identical to the ones without the Casimir energy and given in 



Eqs. ([T9|)-([20]). Therefore, we conclude that the Casimir energy has only sub-leading effect on the radius 
stabilization. 

In [|55ll it is mentioned that the Casimir energy does not affect the radion stabilization due to GW 
mechanism, which is in agreement with our findings. The difference between the two works, is that in 
[|55l the Casimir energy is due to the presence of a conformally massless bulk scalar (since conformal 
invariance is assumed). The effective potential in [|55l behaves as V{a) ~ k^a'^ (1 — da) + 0{a^). 
Stabilization with this potential leads to a tiny mass of the radion m'^^ ~ k'^al which has no effect on the 
GW- stabilization. As with the other stabilization mechanisms based on the Casimir effect, it suffer from 
a fine-tuning problem, [|55l . 



C The Radion Mass from Linearized Einstein Equations 

In this appendix we derive the radion mass in the method of [|50ll , which is extended by [|58l . Here we 
describe trivial extension of it to the flavon case, where there is more than one bulk scalar. Unlike the 
naive ansatz described in Appendix |A], this method takes into account not only the metric's fluctuations 
but also the bulk scalars' fluctuations. In addition, both scalars and metric fluctuations are solutions to the 
linearized Einstein equations. The parametrization used for the radion and the bulk scalars fluctuations is 

Casimir energy ll52l l53l . Since such scalars do not appear in RS phenomenologically viable models and the fermion bulk 
masses are rather model dependent, we neglect their contribution to the Casimir energy. However, bulk gauge fields do appear 
in models were the SM-fields are allowed to propagate in the bulk, therefore we consider their contribution to the Casimir 
energy. 
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ds^ = e'^^^y^-^^^'''y^r]^,dx''dx'' - (1 + G{x, y)fdy'' , (47a) 
(i)i{x,y) = (piiy) + Xi{.x,y) , (47b) 



where A{y) and (piiy) are the background solution, see Eqs. ([9]) and ([T2]). One can add graviton fluctua- 



tions by replacing 77^,^ — )■ r^^;^ + /i^j , were TT denotes transverse traceless. However, as one can check, 
these graviton fluctuations will be decoupled in the linearized Einstein equations (see also [95]). Note 
that in this section we use the notation of [l50l, denoting the bulk potential as V and the brane potentials 
as Auv,iR- The bulk cosmological constant and brane tensions are included in the bulk potential and brane 
potential respectively. From the v components of the Einstein equations one finds the equation 

2d^,d,F-d^d,G = ^. (48) 

Therefore, we set G = 2F from here onwards. This scalar degree of freedom, F, is the radion|^and its 
mass is identified with the lowest eigenvalue of the 4D box operator, □ = rf'^'^d^d^, 

DF = -m^F . (49) 

The linearized Einstein equations, 5Rmn = 5Tm at, where = 1/(2M^). 5_Rm7v is the linearized 
Ricci tensor and 6Tmn denotes linearization of the tensor Tmn = Tmn — \gMNg^^TcD built from the 
energy-momentum tensor Tmn- Their components are 

5R^,u = Vt^uOF + e-^^r]^^{-F" + WA'F' + 6A"F - 2AA'^F) , (50a) 
6R^5 = Sd^F' -QA'd^F, (50b) 
6R55 = 2e^^UF + 4F" - IQA' F' , (50c) 



and 



e 



-2A 



- 3 



mm)x^ - ^vm)F + 5^ {K^m)x^ - ^Km)F) ^{y - y^) 



ritiv^ (51a) 



5T^5 = 0'AX., (51b) 
5T55 = 20:x: + \Vim)x^ + \vm)F + \Y. ^^Umx^ + 2A 5y - y^) , (51c) 

a 

where we use the notation Vi{{(j)\) = - ^ , VjA{(j)\) = S^Y^ 1 1 ^ and similarly for A„; repeated 

* In fact we identify only lowest mass eigenstate as the radion: without the GW mechanism there is only one solution with 
zero mass; with the GW there is a KK tower of solutions for the metric-scalar fluctuations. For details see I.50J. 
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indices are summed over unless otherwise mentioned. The linearized scalar equations of motion are 



- E [^o.^,m)x, + 2A„,({0})F] 5{y - y^) . (52) 

a 

The /i5 equation can be immediately integrated, setting the integration constant to zero 

3F' - 6A'F = K^(t)[xi ■ (53) 

The boundary conditions for the scalars and radion fluctuation are obtained by integrating infinites- 
imally around the delta functions of the /iz/ equation and Xi equation of motion. After use of the back- 
ground jump conditions!^ these are 

[^IL = y i^Umx^ + 2A„({0})F) , (54) 
[x:]L = A,,,({0})x, + 2A,,({0})F. (55) 

Using the background jump conditions the first equation is equivalent to Eq. ( [53] ) and therefore does not 
constrain the solution. 

As in [|50l an eigenvalue equation for F is now formed. This is first done by considering the com- 
bination of e^^(/i, v) + (5, 5) in the bulk, where (A, B) denotes the AB component of the linearized 
Einstein equation. After use of Eq. ( |49l ) one gets 

F" - 2A'F' + 2A"F - SA'^F - e^^m'F = — {2(f)',x'i + AV{{(j)})F) . (56) 

3 

Using Eqs. (|5]) in the bulk further simplifies this to 

F" - 2A'F' - e'^m'F = ^^'.x^ • (57) 



In the case of a single scalar field one can now use Eq. (53 1 to get the eigenvalue equation for F 



F" - 2A'F' - 4:A"F - 2—F' + iA'—F = -m^e^^F . (58) 

(f)' (f)' 

In the case of many scalar fields this elimination is not easily achieved and one generally has to diagonlize 
the entire system (see e.g. Il95l for such attempts). However, for the case considered in section|3]we can 
use the fact that the 0j's have scaled profiles, i.e. 0j = rz/uv,iP(2/), where p{y) is general for all the 
flavons (with z/uv,j = for the fields with vanishing background profiles). In that case we can rewrite 



^ In the notation used here, these are [A'] | = li'^aiW}) ™d [(p'^] | = Xai{{(t>}). 
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Eq. ( [53] ) as 

3 



Therefore, 



^{F'-2A'F)=ruvy,^p\y)x^■ (59) 



^^X'^ = ri^VY.p'x: = -2 - 2A'F) + F" - 2A"F - 2A'F') , (60) 



and using this, one gets Eq. (58 1 with replaced by p. As in [1501 . to find the lightest mode, one solves 



Eq. ( [58] ) perturbatively in r^, using the ansatz 

F{x,y) = e'^''^l + r^f{y))^r{x), m'^ = + r^m\ (61) 

Note that this ansatz is justified by noting that if we neglect A" in the bulk then F = e^^ is a solution of 
zero mass. Denoting 

^(^) = ^ = ^' = % + r2fcG(y) , (62) 

0' P' 



and expanding the Eq. (58 ) to order r^, one gets 



r\y) + 2(A; - Q{y))f\y) - ^k{k - Q{y))G\y) - AkG'\y) + e^^'^m^ = , (63) 



where the boundary conditions given by Eq. ( [55] ). In the case of infinite brane masses the boundary 
conditions are simpler and by use of Eq. ( |53] ) imply 

(F'-2A'F)|^ = 0. (64) 

Expanding this boundary conditions to order yields 

f\y^) = 2kG\y^) . (65) 



The bulk equation Eq. ( [63] ) has two integration constants and the mass parameter, one integration constant 
is an overall normalization factor, while the mass and the other integration constant are found by imposing 
the boundary conditions. We have solved Eq. ( |63| ) and found the radion mass to lowest order in e, r and 
agree exactly with the masses calculated in the effective potential in sections [2| and |3], i.e. Eqs. ( [20| ) 
and([33[). 
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